Kondo effect in carbon nanotube quantum dots with spin-orbit coupling 
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Motivated by recent experimental observation of spin-orbit coupling in carbon nanotube quantum 
dots [F. Kuemmeth et al, Nature (London) 452, 448 (2008)], we investigate in detail its influence 
on the Kondo effect. The spin-orbit coupling intrinsically lifts out the fourfold degeneracy of a 
single electron in the dot, thereby breaking the SU(4) symmetry and splitting the Kondo resonance 
even at zero magnetic field. When the field is applied, the Kondo resonance further splits and 
exhibits fine multipeak structures resulting from the interplay of spin-orbit coupling and Zeeman 
effect. Microscopic cotunneling process for each peak can be uniquely identified. Finally, a purely 
orbital Kondo effect in the two-electron regime is also predicted. 

PACS numbers: 73.23.-b, 73.63.Fg, 72.15.Qm, 71.70.Ej 



Introduction. — The Kondo effect is one of the most fas- 
cinating and extensively studied subjects in condensed 
matter physics Since its first experimental observa- 
tion in semiconductor quantum dots (QDs) in 1998 Q, 
after 10 years of theoretical predictions |3j, various as- 
pects of this effect have been explored in virtue of the 
tunability of relevant parameters in the QDs. The Kondo 
physics was further enriched by fabricating carbon nan- 
otube (CNT) QDs @, H i, 0, IE where additional or- 
bital degree of freedom originating from two electronic 
subbands can play a role of pseudospin. In CNT QDs, 
spin-orbit coupling is widely believed to be weak and 
two sets of spin-degenerate orbits are expected to yield 
a fourfold-degenerate energy spectrum which possesses 
an SU (4) symmetry. Consecutive filling of these orbits 
forms four-electron shells [!, H, 0, 0, S E3] ■ In each shell, 
the SU(4) Kondo effect was observed in the valleys with 
one, two, and three electrons @, 0, H,@]- Theoretically, 
the £{7(4) Kondo effect in CNT QDs has also been ex- 
tensively studied 
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However, recent theories [15| suggest that spin-orbit 
interaction can be significant in CNTs due to their cur- 
vature and cylindrical topology. More recently, transport 
spectroscopy measurements on ultra-clean CNT QDs by 
Kuemmeth et al. [lit ] demonstrate that the spin and or- 
bital motion of electrons are strongly coupled, thereby 
breaking the SU (4) symmetry of electronic states in such 
systems. This motivates us to reconsider the Kondo ef- 
fect in CNT QDs by explicitly taking account of the spin- 
orbit coupling since this symmetry-breaking perturbation 
at the fixed point must break the £[7(4) Kondo effect 
studied previously. An important consequence is that 
even at zero magnetic field the Kondo effect manifests 
as split resonant peaks in the differential conductance. 
At finite fields, these peaks further split into much com- 
plicated subpeaks, reflecting the entangled interplay of 
spin and orbital degrees of freedom. Concerning all mi- 



croscopic cotunneling events involving spin and/or orbit 
flip, these fine multipeak structures can be uniquely iden- 
tified. Moreover, the spin-orbit coupling also determines 
the filling order in the two-electron (2e) ground state [3] , 
producing a purely orbital Kondo effect different from 
that observed by Jarillo-Herrero et al. [6| . 

Model Hamiltonian and QD Green's function. — We 
model a CNT QD coupled to source and drain leads 
by the Anderson Hamiltonian H = Hd + Ht, + fly, 
where H d = J2 m £ ™dl n d m + % J2 m ^ m > n m n m ' , H L = 

J2k, m ,a £ kCt ma C km a, H T = Y,k, m , a V adlnC kma + H.C.. 

Hd models the isolated CNT QD with n m = d m d m . d m 
(d m ) creates (annihilates) an m — {a, A} configuration 
electron in the dot, where a, A = ± are the spin and 
orbital quantum numbers, respectively. U denotes the 
on-site Coulomb repulsion, and e m is the single-particle 
energy. Ht, models the two leads (a = L,R). The tun- 
neling between the dot and the leads is described by Ht 
with spin and orbital conservation since the dot and the 
leads are usually fabricated within a same CNT and thus 
have the same orbital symmetry |(| 0, H, [t| . Assuming 
some orbital mixing, a crossover from SU(4) to SU(2) 
Kondo effect has been investigated [ll|, HH . Here we fo- 
cus on a totally different breaking of the SU(A) Kondo 
effect by the spin-orbit interaction which intrinsically lifts 
out the degeneracy in the single-particle energy. 

The energy spectrum of CNT QDs in low-energy limit 
is Ek — fiVF \k\ with vp the Fermi velocity and k the 
quantized wave vector perpendicular to the CNT axis. 
For a given wave vector fc > 0, k = — Afc (A = ±) de- 
notes two degenerate graphene subbands, corresponding 
to clockwise and anticlockwise classical orbits encircling 
the tube circumference. In curved graphene, the spin- 
orbit interaction is generally classified into three types 
[l5j |. an intrinsic coupling, a Rashba coupling, and a cur- 
vature coupling. However, in most experimentally ac- 
cessible CNT QDs, only the last one is dominant |l5j ]. 
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It can be described by a spin-dependent topological flux 
<J(f>so and shifts k by 6%k = o~4>so I (r4>o) [la |16|, with 
0o being the flux quantum and r the tube radius. A 
magnetic field B introduces an Aharonov-Bohm flux 4>ab 
and further shifts k by 82k = 4>ab I '(x^o), ending up with 
k = — Afco + 5\k + 62k and thus |fc| = fco — XS±k — A^fc- 
Combining \k\ with the linear dispersion, and including 
the spin Zeeman energy, the single-particle energy of the 
state |<7, A) reads e a \ = Ed — erAAso/2 — X^BcosO — oB. 
Here Ed = hvpko is the basic dot level which depends 
on the confining geometry and can be tuned by a gate 
voltage. The second term accounts for the spin-orbit cou- 
pling with Aso — 2frvF(j)so/( r 4 > o)- The last two terms 
represent the orbital and spin Zeeman effects, where 
p = 2norb/(9fJ-B) is the ratio between the orbital mo- 
ment fi or b and the Bohr magneton \ib , the former is typ- 
ically one order of magnitude larger than the latter 17j . 
B = g/isB/2 is the renormalized field, and 9 is the angle 
between the field and the CNT axis. Typical energy spec- 
tra are shown in Fig.[TJ Notably, the spin-orbit coupling 
induces a zero-field splitting Aso and two level crossings 
for the parallel fields at B = ±Bq with Bq = Aso/(2/i). 
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FIG. 1: Ecr\ as a function of (a) parallel and (b) perpendicular 
magnetic fields for Ed = 50meV, p = 10, and Aso = lmeV. 



Electronic transport through the CNT QD is deter- 
mined by the dot retarded Green's function, which is 
G m (e) ee ((d m \dl)) = G m (e)(l + UE m '((n m 'd m \dl))), 
where G° m (e) = [e e m - £° (s)} \ E° ( e ) = £ f 
T = J2 a T a , and T a = irp\V a \ 2 with p being the den- 
sity of states in the leads at the Fermi level. Due to 
the interaction U, G m (e) includes a high-order Green's 
function, whose equation of motion (EOM) produces in 
turn more higher-order ones. We truncate this hierarchy 
by the Lacroix's decoupling procedure [11 03, [I| which 
successfully captures the essential physics of the Kondo 
effect. By this scheme, in the limit of U — » 00, the dot 
Green's function is given by 



G m (e) = 



1 Em' ) ~ Em' 



where the prime in the summation means m' ^ m and 
I 
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[G m > (e')Y 
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with / (e) = (l/r)2 a r a/a(e) where f a (e) represents 
the Fermi distribution function in the leads. Eqs. (JTJ)-((4j) 
can be solved self-consistently. For a QD with an A-fold 
degenerate level, e m — Ed, an SU (N) Kondo temperature 
within EOM is = D exp{7re d /[(A-l)r]} [21], where 
D is the half bandwidth in the leads and Ed is measured 
from the Fermi level. 

Results and discussions. — Experimentally, the mea- 
sured differential conductance dl/dV near zero source- 
drain bias directly reflects the Kondo features in the 
dot density of states Pd{e) = — (1/ir) Im[G m (e)]. 
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Within the Keldysh formalism [22|, the transport cur- 
rent is / = ^^/de \J L {s) - f R (z)}E m MG m (s)}. 
In the numerical results presented below, we use sym- 
metric dot-lead couplings Ti, = Tr and a symmetrically 
applied source-drain bias V. The high-energy cutoff is 
fixed to be D = 1 and T = 0.0LD. 

Figure [2] presents our main results. Even at zero mag- 
netic field, the spin-orbit coupling lifts the degeneracy 
in parallel and anti-parallel spin-orbit configurations of 
a single electron in the dot [see Fig.[T], thereby breaking 
the SUjA) symmetry of the Kondo effect studied previ- 
ously ll|, ll2j, llj, llj] . Instead, the Kondo effect manifests 
as three resonant peaks in the differential conductance, 
which locate at V — 0, ±Aso, respectively. Due to the 
coupling, each peak entangles both spin and orbital de- 
grees of freedom. Specifically the two side peaks arise 
from both spin-flip intraorbital and spin-conserved in- 
terorbital transitions, while the central peak is attributed 
to spin-flip interorbital transitions. A similar zero-field 
three-peak structure has also been predicted in Silicon 
QDs [201 ] but the underlying mechanism is entirely dif- 
ferent. In [2pJ, the two side peaks result from a trivial 
valley (orbital) splitting at zero field and the remaining 
spin degeneracy gives rise to the central peak. 

When the magnetic field is applied, the Zeeman ef- 
fect further removes the remaining degeneracies in the 
energy spectrum and the situation becomes complicated. 
In the field along the tube axis [Fig. [21(a)], owing to the 
interplay of spin-orbit coupling and Zeeman effect, each 
peak at zero field further splits into four subpeaks, ending 
up with rich twelve-peak structures in the Kondo reso- 
nance. At B = Bq, only seven peaks are visible and 
the other peaks merge into the remaining ones. This 
is because of the level degeneracy at this special field 
[see Fig.[lja)]. Though these fine multipeak structures 
appear to be complicated, the inherent physics and the 
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FIG. 2: Fine structures of the Kondo resonance with spin- 
orbit coupling, (a) dl/dV versus V with various parallel 
magnetic fields, B/B = (bottom), 1/3, 2/3, 1, 4/3, 5/3 
(top). The red, green, and blue dashed lines are guides for 
the traces of peak pairs marked by the numbers 1, 2, and 3, 
respectively, (b) dl/dV versus V with various perpendicu- 
lar magnetic fields, B/T^ = (bottom), 0.1, 0.2, 0.4, 0.7, 
1 (top), (c) Schematic representation of three transition pro- 
cesses producing the peak pairs 1, 2, and 3 in (a), respectively. 
In (a) and (b), the curves corresponding to B are offset 
for clarity. The parameters used are Ed = — 10r, /i = 10, and 
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-B-evolution of each peak can be clarified by identifying 
all many-body cotunnelings that add up coherently to 
screen both spin and orbital degrees of freedom. 

In the one-electron regime, the Kondo effect arises 
from the coherent superposition of six transition pro- 
cesses: two spin-flip intraorbital transitions |<r, A) <^> 
— a, A), two spin- flip intcrorbital transitions ex, A) <^> 
| — a, — A), and two spin-conserved interorbital transitions 
| cr, A) | a - , — A). Each transition needs an energy of A, 
(i = 1, ••• ,6) and develops a pair of Kondo peaks at 
V = ±Aj. From energy differences between the initial 
and final states, one readily lists all six transition en- 
ergies as Ai = \e ++ -£-+|, A 2 = |e+_ -e — |, A 3 = 
- £— |, A 4 = |e+_ - e_+|, A 5 = \e ++ - e+_|, and 

A6 = |s i- — e |. Both the spin-orbit coupling and the 

magnetic field are symmetry-breaking perturbations at 
the St/ (4) Kondo fixed point. For Ago = B = and 
hence Aj = 0, the fixed point is reached 11, T^. Thus 
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all the cotunneling processes are elastic and constitute 



FIG. 3: SU(4) Kondo splitting without spin-orbit coupling, 
(a) dl/dV versus V with different parallel magnetic fields, 
B/T^ = (bottom), 0.05, 0.1, 0.2, 0.4, 0.6 (top), (b) 
dl jdV versus V with different perpendicular magnetic fields, 

b / t k ] = (bottom), 0.5, 1, 2, 4, 6 (top). The curves corre- 
sponding to B 7^ are offset for clarity. The parameters used 
are Ed = — 10F, fi = 10, and Aso = 0. 



a highly symmetric St/ (4) Kondo resonance at V — 
(see the curve corresponding to B = in Fig.[3])- On 
the other hand, one can uniquely identify the six tran- 
sition processes from the split Kondo peaks as long as 
all Aj are different from each other and hence the six 
peak pairs (twelve peaks) are well resolved, which re- 
quires Aso / 0, -B ^ 0, and cos# ^ 0. This is exactly 
the case of Fig.^a). As an example, we trace the B- 
evolutions of three peak pairs marked by the numbers 
1, 2, and 3 and unambiguously attribute them to tran- 
sition processes |— ,+) <^ \ — ,— ), | + ,+) \— ,+), and 
| + ,+) | — ,— ), respectively, which are schematically 
shown in Fig.^c). We also note that such an unique 
identification is unavailable in Silicon QDs [20j | where no 
more than nine peaks are visible. 

When the field is applied perpendicularly to the CNT, 
the orbital Zeeman effect is absent and only the spin 
Zeeman effect is involved. As shown in Fig.[2^b), while 
the central peak splits into two subpeaks corresponding 
to the spin-flip intcrorbital transitions, the two side peaks 
split into three subpeaks among which the two split-off 
ones result from the spin-flip intraorbital transitions and 
the rested one corresponds to the interorbital transitions 
without spin-flip. In this case, some peaks are still not 
resolved because A 3 = A 4 and A 5 = A 6 . 

It is useful to comment on the experimental observabil- 
ity of these fine multipeak structures. On the one hand, 
in our EOM scheme the Kondo temperature is underes- 
timated and relatively, the splitting of the Kondo peaks 
is more evident. On the other hand, the decoherence 
[23| neglected in the present study has an effect to smear 
out the split peaks, especially those with high energies. 
These two features will make the observation a bit diffi- 
cult. However, in ultra-clean CNT QDs and using highly 
resolved spectroscopy measurements, it is still possible to 
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observe part or all of these fine structures. 

For comparison, we plot in Fig.[3]the multiple splitting 
of the 5/7(4) Kondo resonance by neglecting the spin- 
orbit coupling. On increasing the field, the SU(4) Kondo 
peak at the zero bias splits in a simple way following the 
Zeeman effect, which produces characteristic structures 
in agreement with the results in the literature @, 
1 2( | but quite different from that discussed above. 
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orbital Kondo effect develops in the 2e ground state due 
to the particular multielectron filling order. Our results 
indicate that the spin-orbit coupling significantly changes 
the low-energy Kondo physics in CNT QDs. 

Support from the NSFC (10575119), the Major State 
Basic Research Developing Program (2007CB815004), 
and the Program for NCET of China is acknowledged. 




FIG. 4: Purely orbital Kondo effect in the 2e regime, (a) 
dl/dV versus V with various parallel magnetic fields, (B — 
B )/T ( k ] = -0.6 (bottom), -0.4, -0.2, 0, 0.2, 0.4, 0.6 (top). 
The curves are offset for clarity and the parameters used are 
e' d — — 3F and fi — 10. (b) Schematic diagram of the inherent 
many-body cotunncling. 

The spin-orbit coupling further determines the filling 
order in the 2e ground state [lj|, which exactly fol- 
lows the first excited state of the single-particle spectrum 
(Fig. [I}. At the parallel fields near the degenerate point 
B , while the first electron always occupies the state 
| + , +), the second electron can fluctuate between | — , — } 
and | — ,+}, giving rise to a purely orbital Kondo effect 
as its spin a = — is fixed. To describe this Kondo effect, 
m appearing in all our previous formulae should be re- 
placed by { — , A} and the single-particle energy becomes 
e'_ x = e-\ + U 2e = e' d +(l- \li)(B-B ), where U 2e is the 
2e Coulomb interaction and e' d — Ed + Aso/(2/^) + Ui e ,- 
Fig-SJa) presents the resulting dl/dV as a function of V. 
For B = B , there is a pronounced zero-bias peak which 
represents an SU (2) orbital Kondo effect and provides a 
spin-polarized conducting channel. The peak splits due 
to the field applied away from the degenerate point. This 
orbital Kondo effect, as schematically shown in Fig.[4jb), 
comes from the same shell and dwells in the 2e valley, 
therefore being distinct from the one recently observed 
in CNT QDs without spin-orbit coupling 

Conclusion. — We have studied the Kondo effect in a 
CNT QD with spin-orbit coupling. It is shown that the 
Kondo effect manifests as rich fine multipeak structures 
in the differential conductance when a magnetic field is 
applied. These fine structures are quite different from the 
SU(A) Kondo effect studied previously and might be ob- 
servable in future experiments. In such a system, a purely 
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